The relative importance of electron-electron interactions compared to disorder in the 

two-dimensional "metallic" state 



O 

o 
< 

in 



CD 



A. Lewalle ,^ M. Pepper/ C. J. B. Ford,^ E. H. Hwang, ^ S. Das Sarina,^ D. J. Paul,i and G. Redmond^ 
^ Cavendish Laboratory, Madingley Road, Cambridge CB3 OHE, United Kingdom 
^ Department of Physics, University of Maryland, College Park, Maryland 20742-4111 
^ Nanotechnology Group, National Microelectronics Research Centre, Lee Maltings, Prospect Row, Cork, Lreland 

(February 1, 2008) 

The effect of substrate bias and surface gate voltage on the low temperature resistivity of a 
Si-MOSFET is studied for electron concentrations where the resistivity increases with increasing 
temperature. This technique offers two degrees of freedom for controlling the electron concentration 
and the device mobility, thereby providing a means to evaluate the relative importance of electron- 
electron interactions and disorder in this so-called "metallic" regime. For temperatures well below 
the Fermi temperature, the data obey a scaling law where the disorder parameter (kpl), and not the 
concentration, appears explicitly. This suggests that interactions, although present, do not alter the 
Fermi-liquid properties of the system fundamentally. Furthermore, this experimental observation is 
reproduced in results of calculations based on temperature-dependent screening, in the context of 
Drude-Boltzmann theory. 

PACS numbers; 71.30.-t-h, 73.40. Qv 
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The term "metallic behaviour" , for a two-dimensional 
(2D) system with concentration n greater than some crit- 
ical value He, has come to designate a drop in the resis- 
tivity p as the temperature T is decreased, as opposed to 
"insulating behaviour" (for n < Uc) for which p increases 
exponentially with decreasing temperature. In experi- 
ments performed in systems with low disorder and small 
values of rtc p|-p^ , the change in p in the metallic regime 
is considerably greater than was seen in earlier work at 
higher densities and in more disordered systems |pT| , |l2| . 
There is still no consensus on whether this change in p 
is caused by Coulomb interactions or disorder. The low 
values of ric (oc r~^, where rg is the ratio of the Coulomb 
interaction to Fermi energy), the scaling behaviour of 
transport properties around ric and the apparent satura- 
tion of p(T) (with no sign of an upturn), seen by some 
authors P,p|,p|, p^|jT^ [l5t as T 0, have led to sugges- 
tions that the system is a true metal (in the sense that 
the charge carriers are delocalized at T = 0). The sug- 
gestion is that this novel phase results from strong elec- 
tronic correlations, and that ric marks a quantum critical 
point and a metal- insulator transition (MIT). This con- 
tradicts the well-established principle that a 2D system 
of non-interacting carriers is insulating at T = p6[ , but 
there is no definitive corresponding prediction when in- 
teractions are present. Other authors, however, observe a 
negative magnetoresistance at low magnetic fields B and 
weak-localization corrections to the resistivity, which are 
consistent with Fermi-liquid behaviour [p7|-pO{. Many 
hypotheses have been put forward to explain the metal- 
lic behaviour, some invoking exotic interaction effects 
and others advocating a more traditional frame- 
work, suggesting that the "metallic" behaviour is only a 
finite-temperature effect that is overwhelmed by localiza- 



tion at sufficiently low temperatures ]24|-|29||. To a large 
extent, the debate hinges on the question of whether the 
relatively large values of Tg (typically > 10), correspond- 
ing to Uc, warrant fundamentally new physics. 

In motivating our work, we note that as n (rg) de- 
creases (increases), it is not only the electron-electron 
interaction that is increasing in the 2D system, but also 
the effective disorder felt by the carriers. This is be- 
cause the dominant disorder in semiconductor 2D sys- 
tems arises from Coulomb scattering by charged impuri- 
ties, which increases monotonically with decreasing den- 
sity as the 2D carrier system becomes less effective in 
screening the Coulomb interaction between the carriers 
and the charged impurities. Since, in the absence of umk- 
lapp processes (not applicable in these systems) , electron- 
electron interactions typically do not affect ohmic trans- 
port, there is good reason to believe that much of the ob- 
served "metallic" behavior may be arising from the weak- 
ening of screening in the disorder potential, rather than 
from the increasing electron-electron interaction. By 
carefully studying temperature-dependent 2D transport 
in high quality Si-MOSFETs, where disorder strength 
and carrier density are controlled independently using 
a substrate bias, we hope to shed light in this Letter 
on this question of the relative importance of electron- 
electron interactions and the disorder potential in low 
density "metallic" transport. 

A negative substrate bias (Kiub) steepens the trian- 
gular confining potential normal to the surface, pushing 
the electron wave-function closer to the Si-Si02 inter- 
face. This can either increase or decrease the mobil- 
ity 1^,^,^. The surface-gate voltage is used to tune 
n. Without assuming either a specific mechanism for 
the increase in p with temperature or the existence of a 
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quantum phase transition at Uc, this technique probes 
the density-mobihty phase space in a way that is im- 
possible with singly-gated devices, and it allows a study 
of the metal in terms of Tg and the disorder parame- 
ter fcp/, where fcp is the Fermi wave- vector and / is the 
momentum relaxation length. The value of kpl can be 
calculated directly from the relation p = {h/2e'^)/ (kpl), 
whilst n is measured from the Hall effect at magnetic 
fields B < 0.5 T. 
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FIG. 1. Threshold voltage and peak mobility versus ap- 
plied substrate bias for T=1.4 K. 



Measurements were performed on an n-type Si- 
MOSFET inversion layer with a 200 nm-thick oxide layer, 
fabricated by NMRC, with a peak mobihty (at T = 1.4 K 
and with V^ub = 0) of 1.9 m^V^^s"^ corresponding to 
n = 4 X 10^^ m^^ on 100 flcm substrate. Devices were 
Hall bars of dimension 1000 x 100 pm^. A "^He cryo- 
stat was used, covering a range 1.4 K < T < 70 K. 
Resistivities were measured using a standard a. c. four- 
terminal technique, with a constant source-drain current 
of 100 nArms and a frequency of 19 Hz, after ensuring 
that the magnitude of the current did not induce signif- 
icant electron heating above the substrate temperature 
and that the contacts were ohmic. A substrate bias Vsuh 
was applied at room temperature before the sample was 
lowered into the cryostat and cooled slowly. 

One consequence of the (negative) substrate bias is to 
increase the threshold voltage, which results in the in- 
crease of the confining electric field, for a given n. The 
peak mobility (to be considered here as a rough measure 
of the disorder) decreases as Vsuh is made more negative, 
as shown in Fig. ^ consistent with the enhancement of 
scattering from interface roughness and the charged im- 
purities located near the Si02 interface. The enhance- 
ment of spin-orbit scattering is another possible conse- 
quence. Negative magnetoresistance is observed in a 
weak magnetic field perpendicular to the system, how- 
ever, in agreement with the prediction for the quenching 
of weak localization. No sign of the positive magnetore- 
sistance associated with spin-orbit scattering is seen , 
indicating that this is not a strong effect. 




FIG. 2. Phase decoherence time versus T'^ for a concen- 
trations 1.56, 3.48, and 4.42 x 10^^ m"^ at V^ub = 0. 

The magnetoresistance demonstrates that weak lo- 
calization, which results from quantum interference, is 
present in spite of the metallic behavior. Fig. ^ shows 
that the phase decoherence time r^, obtained by fitting 
the negative magnetoresistance, is proportional to T~^, 
in agreement with Fermi liquid theory ||3^ , |33|] . At suffi- 
ciently low temperatures, it is likely that this localization 
will dominate. 
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FIG. 3. p(r) over a range of densities for three substrate 
biases —20, —40, and —60 V (solid curves). Dotted curves 
show p(r) for Vsub = 0. Ranges of densities are indicated 
by bold numbers (xlO^^ m"^) for Kiub < 0. Numbers in 
parentheses give the density range of the Kub = curves. 
The critical concentration ric ranges from 0.99 x 10^^ m~^ (for 
= 0) to 1.41 X 10^^ m-2 (for Kub = - 60 V). Fermi 
temperature Tp [K] = 7.25 n [10^^ m*^]. 

The solid curves in Fig. ^ show p(r) for different values 
of Vsub and for a range of densities spanning the "metal" - 
insulator transition. They are superposed on dotted 
curves, which correspond to V^ub — 0. In each figure, the 
transition from strong localization {dp/ dT < 0) to metal- 
lic behavior occurs at the same value of the resistivity (or 
kpl), pc « h/e"^, whereas ric varies from 0.99 x 10^^ m~^ 
(at y^ub = 0) to 1.31 X 10^5 m-2 (at V;ub = -60 V). 
This critical resistivity corresponds to kpl = 0.5, a value 
below which the Fermi wavelength is poorly defined. It 
is worth emphasizing that it is kpl and not which is 
the critical parameter. For densities greater than those 
shown in Fig. | {i.e. n > 2 x 10^^ m'^), p{T) IS non- 
monotonic, bending downwards in an insulator-like fash- 
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ion as T is increased beyond Tp, the Fermi temperature. 
It also shows signs of saturation as T approaches zero, 
possibly due to collision broadening or the finite Dingle 
temperature (discussed further below) . The intermediate 
temperature range shows an approximately linear depen- 
dence. The gradients of the curves with T^ub < differ 
from those of the Vsnh = reference curves, showing 
that the substrate bias does not only amount to a rela- 
belling of the p{T) curves. 

With the hope of finding some underlying universal- 
ity, we follow the work of A. R. Hamilton on GaAs [ p4[ , 
and consider theoretical studies, done in the context of 
Drude-Boltzmann theory, which take into account the 
strong temperature dependence of screening |26| , |35| , p6| . 
These calculations predict metallic behaviour resulting 
from the anomaly in the polarizability function n(q) at 
a value of the scattering wave vector q — 2kp which 
dominates the mobility. To lowest order , 

p(r;n) = p(T = 0;n)[l + C(n)T/TF] (1) 

vahd in a range T/Tp < 1 < T/Td, where Td is the ef- 
fective Dingle temperature and C is a function of n and 
of the dominant scattering mechanism. 
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FIG. 4. (a) Experimental p{T /Ty) for a range of densities 
at a constant substrate bias, Vsub = 0. Tp ranges from 11.3 K 
to 38.0 K. The effective values of T/Tp displayed include the 
correction for collisional broadening (see Eq. ^ . (b) The same 
data with the vertical axis scaled by a factor po determined 
empirically for each curve. 

Following Eq. ^, the temperature axis is normalized 
as T/Tp, where T-p is the Fermi temperature calcu- 
lated for each curve individually, as shown in Fig. ^(a), 
for Vsnh = 0. The Dingle temperature was esti- 
mated from Tb = where Tq is the quantum 
lifetime, measured from the amplitude of Shubnikov-de 
Haas oscillations [p^. For a range of concentrations 
n = 4.6 - 5.4 X lO^'^ m-2 (Tq ~ 35 K), we obtain 
Td ~ 1-1 K. At higher T, the curves are linear with 
a gradient that increases with decreasing n. It is remark- 
able that another scaling factor po j applied to the vertical 
axis and determined by eye for each curve, succeeds in 



collapsing the linear portion of all of the curves onto the 
bottom (high-density) curve (Fig. |^b), at least for a range 
of T well away from Tp. Indeed, this scaling factor must 
equal the ratios of both the gradients of the lines and their 
intercepts: 

p(T;n)/po(n)-/(r/rF); T/Tf<0.5. (2) 

This simple transformation holds well for n > 3 x 
IQis m-2 (for which p{T = 1.4 K) < 1.4 m/U) but 
fails at densities closer to the MIT, i.e. it is impossible 
to match both the gradients and the intercepts with a 
single scaling factor. A much better scaling of the curves 
is recovered by including a correction: 

(r/rF)effoctivc = [(T/Tf)2 + A/{kplfY'\ (3) 

The value of the coefficient A was set to 1.0, but the over- 
all quality of the fit was not sensitive to its exact value. 
The correction term 1/k-pl may be interpreted as repre- 
senting the broadening of the Fermi circle due to the finite 
scattering length (collisional broadening) |^^. Alterna- 
tively, it may be viewed as a lowering of the effective den- 
sity (or Tp) resulting from carrier freeze-out, a possibility 
envisaged by some authors [^,^ . With this correction, 
Eq. ^ is obeyed down to n = 1.56 x 10^^ m~^ (for which 
p{T — 1.4 K) = 8 kil/n), i.e. considerably closer to the 
MIT. At even lower values of n, the correction term is 
much more significant compared with T /T-p, so that the 
applicability of Eq. ^ becomes questionable. Neverthe- 
less, the success of Eq. ^ at slightly larger concentrations 
suggests that the low-n regime close to the MIT is where 
the interaction between screening and disorder dominates 
and cannot be treated by a simple perturbative analy- 
sis. For the higher-density data (n > 3 x 10^^ ni~^), 
the broadening correction to the temperature requires 
small modifications to the empirical po values in order 
to achieve a satisfactory collapse, but the quality of the 
collapse is not significantly altered. 
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FIG. 5. (a), (c) p{T/Tf) for fixed densities 2.30 x 10^^ m"^ 
(Tp = 16.7 K) and 4.45 x 10^^ m'^ (Tp = 32.3 K) for sub- 
strate biases 0, —20, —40, and —60 V. (b) and (d) show the 
same data after scaling by po(V^ub). 

Fig. ^ shows the same procedure, applied to a set of 
curves corresponding to different values of V^ub but with 
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the same density (2.30 and 4.45 x 10^^ m~^). Again, 
Eq. |2| is obeyed, but with po = po(Vsub)- The similarity 
in the forms of Eqs. |^ and ^ is notable, but the scaling 
factors Po of Eq. || may be equated with the p{T = 0) of 
Eq. |l| only in so far as (a) p{T] n) remains linear outside 
the strict range T <$:Tf, (b) the saturation observed at 
low temperatures may be discarded as due to the finite 
Td, and (c) C is a sufficiently weak function of n. Ex- 
perimentally, a constant C « 8 is observed, considerably 
greater than the expected theoretical values 2.0 - 3.0 p6[ . 
In spite of this discrepancy, the experimental scaling be- 
haviour suggests that the value of p (or kpl) at finite T 
is determined solely by its value at low T, with the con- 
centration (and therefore Tg) appearing only implicitly in 
Tp. Combining Eq. ||with the identity p = {h / ie^) / {k^l) 
we obtain 



kFl{T ^ 0) ^ /(T/Tp) 
/(O) 



(4) 
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FIG. 6. (a) Numerical calculation of p{T /Tf) due to tem- 
perature-dependent screening, for constant 
n = 2.3 X lO^'^ m"^ (Tp = 16.7 K) and Nu = 1.5, 2.5, 4.5, 8.0, 
and 15.0 xlO^'^ m'^ withTo/rp = 0,0.06,0.12,0.18, and 0.24 
respectively, (b) The same data, with p scaled by p{T = 0). 

Fig. ^ shows numerical results of a numerical calcu- 
lation of the temperature-dependent screening for ex- 
perimental conditions analogousto those in Fig. |^(a,b). 
It includes the temperature dependence of Ti{q) at fi- 
nite temperature, the finite extent of the wavefunction 
in the direction normal to the interface (using a Fang- 
Howard distribution |38|), and a variable effective dis- 
order, parametrized by the the collision-induced Dingle 
temperature Td- (Here, Td measures the strength of 
impurity scattering and is, by definition, temperature- 
independent.) Carrier freeze-out [ p7| was not assumed. 
The penetration depth of the wave function below the 
interface is determined by the depletion concentration, 
iVD, which increases as T4ub is made more negative. The 
range of N^i indicated in Fig. ^ is consistent with the 
range of Vsuh (0 to —60 V) used in the experiment Q. 
By setting suitable values of Td, the calculation repro- 



duces the experimental data with reasonable accuracy: 
(1) p(T) increases in a metallic fashion, with a weaker 
temperature dependence appearing (at low T) as T/Td 
decreases and (at high T) as T/Tp increases; (2) The 
scaling law, Eq. ^, is obeyed (with po = po(A^d)), al- 
though the ability to scale the theoretical curves as was 
done with the experimental data probably relies on the 
correct combination of A^d and Td. We emphasize, how- 
ever, that the theoretical results presented here are not 
intended to give a fully quantitative description of the 
experimental results. Rather, we demonstrate that, to 
a large extent, temperature-dependent screening effects, 
in a conventional transport theory, can account for the 
main features of p(T) in the so-called metallic regime in 
the absence of a magnetic field. 

In summary we have demonstrated that a simple em- 
prical scaling law, where the concentration does not ap- 
pear explicitly, is applicable on the metallic side of the 
MIT when either the concentration or the mobility is 
varied. (A similar effect has been observed in GaAs 
heterostructures The quality of the scaling is im- 

proved, and is obeyed down to lower concentrations, by 
taking into account collision-broadening effects which en- 
hance the effective temperature. This suggests that elec- 
tronic interactions are not the most crucial element in 
the metallic phenomenon, apart from screening the scat- 
tering potential, and that the "metallic" state is basically 
a Fermi liquid. In fact, we observe a striking similarity 
with the prediction of Drude-Boltzmann transport the- 
ory in the presence of temperature-dependent screening 
(see also Ref. [|l9|). Although quantitative agreement be- 
tween this theory and our experimental results is only ap- 
proximate, it is clear that many features of the "metallic" 
behavior can be understood in terms of screening and dis- 
order, without resorting to electron-electron interaction 
or quantum interference effects. 

This work was funded by the UK EPSRC. 
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